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The MC@NLO method as implemented in the Sherpa MC generator is presented using the 
production of W-bosons in conjunction with up to three jets as an example. Corresponding 
results computed at next-to leading order in QCD and including parton shower corrections 
are compared to recent experimental data from the Large Hadron Collider. 



1 Introduction 

To make the LHC a discovery machine we have to acknowledge the fact that it is a QCD 
machine. Many signals suffer from large backgrounds largely due to QCD multi-jet production 
which have to be under good theoretical control to interpret the measurements. 

There are mainly two approaches to include higher-order QCD corrections in theoretical 
calculations of scattering matrix elements (ME): 

Fixed-order ME calculations put an emphasis on the exact evaluation of a hnite number 
of terms in the perturbation series. Apart from being exact to the given order this also has the 
benefit of including all interference terms from different diagrams and working with a number 
of colours Nc — 3. Last, but not least, it becomes possible to include also the exact finite part 
of virtual contributions in a fixed-order calculation. 

Their disadvantages appear when an observable becomes sensitive to logarithmically en- 
hanced regions. It is not sufficient to truncate the perturbation series at a fixed order anymore, 
if the perturbative nature of the coupling constant a s is countered by large logarithms which 
appear when partons become soft or collinear to each other. This problem is solved in the par- 
ton shower approach (PS), where the logarithmically enhanced contributions are resummcd 
to all orders, albeit only in an approximation valid in the collinear limit of the matrix clement 
and in the large N c limit. This allows to generate events with partons at the hadronisation 
scale and thus enables exclusive hadron- level event generation. 

It is thus a worthwhile goal to combine the two approaches and keep the advantages: Include 
the virtual contributions and first hard emission from the exact ncxt-to-lcading order matrix 
element, and add further parton evolution using a parton shower approach. 
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2 Recap: Resummation and NLO calculations 



The basic property of QCD allowing a parton shower resummation is the universal factorisation 
of real emission matrix elements in the collincar limit: 
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With the approximation that multiple emissions happen independently of each other (thus 
yielding Poisson statistics) the corresponding branching probability can be exponentiated to 
give the total no-branching probability down to an evolution scale t: 
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To understand the implications of the no-branching probability A, let us look at the expec- 
tation value of an observable O taking into account up to one emission from the parton shower 
on top of a Born-level event: 
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The "unresolved" contribution corresponds to the event generation case where no emission 
above the parton shower cut-off scale to has been generated and is thus proportional to the 
no-branching probability A (to). The "resolved" contribution on the other hand represents the 
integration over events which had an emission with evolution scale t > to- 

As a reminder and to fix some notation, the contributions of an NLO calculation for the 
expectation value of O are briefly recalled: 
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Here, the Born (£>), virtual (V) and real emission (TV) pieces are denoted accordingly. Since 
V and 1Z are oppositely divergent in separate phase space integrations, one needs to employ 
a subtraction procedure: The real subtraction terms T> are linked to their integrated form X 
by a phase space integration over the 1-emission phase space and can be calculated e.g. in the 
scheme of p]. 



3 Combining NLO and PS 

Applying PS resummation on top of such NLO events is not straightforward: The observable 
has to be evaluated at different kinematics in the 1Z and V terms. But if that is taken into 



2 



DIS 2012 



account, and they are thus showered separately, one generates an additional term at 0(a s ) [2] 
which is not present in the NLO calculation ("double counting"). To counter that term a 
solution was proposed in [2] introducing an additional set of "modified" subtraction terms 
X>( A ). When generating events according to that modified NLO cross section, they will have 
either kinematics (resolved, non-singular term) and are kept as they are or $s kinematics. 
In the latter case, they are processed through a one- step PS with A( A \ i.e. using the modified 
subtraction terms as PS kernels, either generating an emission (resolved, singular) or no emission 
(unresolved, singular) above t . The result of this procedure, 
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with £ (A) ($s) defined as 



z? (A) (<i> s ) = S($ s )+V($ B )+^xf ) ( ( i> B )+^ 



{«} fi=q,g 



can be shown to reproduce (0)( NL0 ) to 0(a s ). 

This procedure still leaves the freedom of choosing A^ A ). The original approach [2] uses the 



parton shower splitting kernels as additional subtraction terms, V 



(A) 



V 



(PS) 



This option 



has the advantage that the exponentiation in the "resolved, singular" contribution is naturally 
bounded by the factorisation scale /Up. Problems appear though with soft divergences in the 
"resolved, non-singular" integration, since the parton shower splitting kernels do not cover soft 
divergences. 

An alternative approach was suggested in [3] and implemented in SHERPA |4j , where the full 
Catani-Seymour dipoles are used pj^ — > T>^' . With this, simplifies significantly, but at a 
cost: can become negative, resulting in A > I. This is generated in SHERPA by a weighted 
Nc = 3 one-step PS based on the subtraction terms V^ S K With this approach, exact NLO 
accuracy also for sub-leading colour configurations is achieved. The phase space boundary for 
the exponentiation though has to be imposed "artificially" by cuts in the dipole phase space. 



4 Results 

Results for W + n-jet production at the LHC are presented here in comparison to ATLAS 
data H] . Events are simulated using SHERPA's Mc@NLO for W + 0, W + l, W + 2 and VF + 3-jet 
production with virtual corrections from BLACKHAT [Sj including a leading-colour approxima- 
tion for the W + 3-jet virtual ME. For n > events are generated requiring fcy jets with 
p± > 10 GeV, and the exponentiation region was restricted using an a = 0.01-cut in the dipole 
terms [7]- The comparison comprises of three levels of event simulation: "NLO" as fixed-order 
calculation, "MC@NLO lem" as MC@NLO including the hardest emission only and "MC@NLO 
PL" as Mc@NLO including the full PS. All distributions are simulated at NLO accuracy and 
the theoretical predictions describe data very well. 



DIS 2012 



3 



Jet transverse momenta Rapidity Distance of Leading Jets 




o . I I I I I I I I I I I I I I I I I I I I I I I I L 




50 100 150 200 250 300 og t u5 2 z5 3 

P ± [GeV] A^(First Jet, Second jet) 



Fig. 1: Transverse momentum of the first, sec- Fig. 2: Angular correlations of the two leading 
ond and third jet (from top to bottom). jets in W 1 * 1 -)- > 2 jet production. 



5 Conclusions and outlook 

The concept of NLO+PS matching was presented in a common formalism, and Mc@NLO was 
developed as a special case. It was emphasised that an exact treatment of sub-leading colour 
configurations can be achieved by exponentiating Catani-Seymour subtraction terms. The first 
NLO+PS predictions for VF+3 jets were presented and showed good agreement with exper- 
imental data from ATLAS. With this method at hand, it becomes now possible to generate 
exclusive hadron-level event samples at NLO accuracy. 

In the future, an improved functional form of the phase-space cut a will allow for a better 
limitation of the exponentiation region. Also, by merging NLO+PS with higher-multiplicity 
tree-level MEs one can already provide an improved description of multi-jet final states in the 
MENLOPS approach [5]. Ultimately, it remains a goal to achieve the merging of NLO+PS 
predictions at different multiplicities into one inclusive sample. 
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